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Abstract 

We study the branching ratios of the lepton flavor violating processes [i — > e^y, r — > ej 
and t — > in the split fermion scenario, in the framework of the two Higgs doublet model. 
We observe that the branching ratios are relatively more sensitive to the compactification 
scale and the Gaussian widths of the leptons in the extra dimensions, for two extra 
dimensions and especially for the r — > ^7 decay. 
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1 Introduction 



The radiative lepton flavor violating (LFV) decays exist at least at one loop level and they are 
rich from the theoretical point of view since the measurable quantities of these decays contain 
number of free parameters of the model used. In framework of the standard model (SM), these 
LFV decays are forbidden and, by introducing the neutrino mixing with the non zero neutrino 
masses, the LFV decays are allowed. However, their branching ratios (BRs) are tiny and much 
below the experimental limits due to the smallness of neutrino masses. This forces one to go 
beyond the SM and LFV decays are among the strong candidates for testing the possible new 
physics effects. 

There is an extensive experimental and theoretical work done on the radiative LFV decays 
in the literature. The current limits for the branching ratios (BRs) of fi — » ej, r — > ej and 
r — > fi'y decays are 1.2 x 10~ n [lj, 3.9 x 10~ 7 [2] and 1.1 x 10~ 6 [3j, respectively. The theoretical 
analysis of these decays has been performed in various studies [I]- [15]. [I]- [10] were devoted to 
the analysis in the supersymmetric models. In [TTj [T2l [T3| ITS] and [15], they were examined 
in the framework of the two Higgs doublet model (2HDM) and in a model independent way 
respectively. With the extension of the Higgs sector and the assumption that the flavor changing 
neutral currents (FCNC's) at tree level are permitted, the BRs of the radiative decays under 
consideration can be enhanced theoretically up to the experimental limits. 2HDM with tree 
level FCNC currents is one of the candidate, and, in this model, the radiative LFV decays are 
induced by the internal new neutral Higgs bosons h° and A . 

This work is devoted to the LFV processes [i —>■ ej, r — > e'j and r — > /ry in the 2HDM, 
with the inclusion of one (two) extra spatial dimension. Here, we choose that the hierarchy 
of fermion masses is coming from the overlap of the fermion Gaussian profiles in the extra 
dimensions, so called the split fermion scenario [lEJdT]. In this case the fermions are assumed 
to locate at different points in the extra dimensions with the exponentially small overlaps of 
their wavef unctions. There are various works done on this scenario in the literature [18]- [35]. 
The explicit positions of left and right handed components of fermions in the extra dimensions 
have been predicted in [18] . Using the leptonic W decays and the lepton violating processes, 
the restrictions on the split fermions in the extra dimensions have been obtained in [19]. The 
CP violation in the quark sector has been studied in [2U] and to find stringent bounds on the 
size of the compactification scale 1/R, the physics of kaon, neutron and B/D mesons has been 
analyzed in [21] . [22] is devoted to the rare processes in the split fermion scenario and [231 121] 
is related to the shapes and overlaps of the fermion wave functions in the split fermion model. 
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In [25] the electric dipole moments of charged leptons have been predicted, respecting this 
scenario. 

In the present work, we study the BRs of the radiative LFV decays by considering that the 
leptons have Gaussian profiles in the extra dimensions. First, we consider the BRs in the case 
of a single extra dimension. In the following, we make the same analysis when the number 
of extra dimensions is two and the charged leptons are restricted to the fifth extra dimension, 
with non-zero Gaussian profiles. Finally, we estimate the effects of the extra dimensions if 
the non-zero Gaussian profiles exist in both extra dimensions. We observe that the BRs are 
relatively more sensitive to the compactification scale and the Gaussian widths of the leptons 
in the extra dimensions, for two extra dimensions and especially for the r — ► /ry decay. 

The paper is organized as follows: In Section 2, we present the BRs of the radiative LFV 
decays in the split fermion scenario, in the 2HDM. Section 3 is devoted to discussion and our 
conclusions. 

2 The radiative lepton flavor violating decays in the split 
fermion scenario in the two Higgs doublet model 

The radiative LFV decays are worthwhile to study since they exist at least in the loop level 
and they are rich theoretically. The tiny numerical values of the BRs of these decays in the SM 
forces one to go beyond and the version of the 2HDM, permitting the existence of the FCNCs 
and the LFV interactions at tree level, is the simplest candidate. The new Yukawa couplings, 
which are complex in general, play the main role in the calculation of the physical quantities 
related to these decays. In addition to the extension of the Higgs sector, the inclusion of the 
spatial extra dimensions brings additional contributions. Here, we take the effects of extra 
dimensions into account and we follow the idea that the hierarchy of lepton masses are coming 
from the lepton Gaussian profiles in the extra dimensions. The Yukawa Lagrangian responsible 
for these interactions in a single extra dimension, respecting the split fermion scenario, reads 

£ Y = Zi ij liL<t>2E jR + h.C. , (1) 

where L and R denote chiral projections L(R) = 1/2(1 =p 75), <p 2 is the new scalar doublet and 
£5^ - are the flavor violating complex Yukawa couplings in five dimensions. Here, Ul (Ejn), with 
family indices i, j, are the zero modqj lepton doublets (singlets) with Gaussian profiles in the 



1 Notice that we take only the zero mode lepton fields in our calculations. 
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extra dimension y and they read 



E jR 



it i 



N e -(y~yjn) 2 / 2 ° 2 e 



jR, 



(2) 



with the normalization factor iV = x/ i 1/2 ■ Ul (EjR) are the lepton doublets (singlets) in four 
dimensions and the parameter a represents the Gaussian width of the leptons, satisfying the 
property a « R, where R is the compactification radius. In eq. (j2J), the parameters yu and 
Dm are the fixed positions of the peaks of left and right handed parts of i th lepton in the fifth 
dimension and they are obtained by taking the observed lepton masses into account [IE]. The 
underlying idea is that the mass hierarchy of leptons are coming from the relative positions 
of the Gaussian peaks of the wave functions located in the extra dimension [T6J [T7J [18]. One 
possible set of locations for the lepton fields in the fifth dimension read (see [TS] for details) 
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We choose the Higgs doublets 0i and 2 as 
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with the vacuum expectation values, 

< 0i >= 
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(5) 



and collect SM (new) particles in the first (second) doublet. Notice that Hi and H2 are the 
mass eigenstates h° and A respectively since no mixing occurs between two CP-even neutral 
bosons H° and h° at tree level in our case. The LFV interaction at tree level is carried by the 
the new Higgs field 02 and it is accessible to extra dimension. Following the compactification 
on the orbifold S 1 /Z 2 , it reads 

1 



2 O) (x) + v / 2£4 n) W <x*{ny/R) 



(6) 



where 02°' ) (^) (02^ i x )) is the Higgs doublet in the four dimensions (the KK modes) including 
the charged Higgs boson H + (H^ + ), the neutral CP even-odd Higgs bosons h°- A (/i 0(n )- 



rn 



H± 



+ mi , and the neutral 



A°( n >). The non-zero n th KK mode of the charged Higgs mass is 
CP even (odd) Higgs mass is y^fo + tr^, ( \J m \o + m n )> with the n'th level KK particle mass 
m n = n/R. 
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The radiative decays under consideration exist at least in the one loop level with the help 
of the intermediate neutral Higgs bosons h , A and their KK modes (see Fig. [1]). The 
lepton-lepton-S" vertex factors V^fRn^ in the vertices f iL ^ S^ n '(x) cos(ny/R) fj R (L), with 
S = h°,A° and the right (left) handed i th flavor lepton fields fj R (L) in five dimensions (see eq. 
(j2J)), are obtained by the integration over the fifth dimension and they read 

V£ R(RL)ij = e~ n2 ° 2 / iR2 e-^W-y^ 2 ^ 2 cos [ w (^W+^w) ] . (7) 

In the case of n = 0, the factor becomes Vlrirl)^ = e~ ( ~ yiL< - R ^~ yiR< - L ^ 2 ^ 4:a2 and we define the 
Yukawa couplings in four dimensions as 

e| ((£f) f ) = vE R{RL)ij ^ ((^/V^r . (8) 

Since the LFV processes , fi — ■> ej, r — > and r — > /ry exist at loop level, there appear 
the Logarithmic divergences in the calculations and we eliminate them by using the on-shell 
renormalization scheme Q Taking only r lepton for the internal line the decay width T, 
including a single extra dimension, reads as 

r(/i-/ 27 )=c 1 (|A 1 | 2 + |A 2 | 2 ) , (9) 

where 

Ai = Qt— — j ( 6 m r ln*rh ^N*hr ( F ( z h°) - F{z A o) + 2 ^ c' n (A, t) c n (/ 2 , r) 

46m r \ rt=l 

x (F(z„, ft o) - F^o))) + m A e^ r/2 £$, t/i (G(zfco) + G(z A o) 

oo \ 

+ 2 £ e-™ 2CT2 / 2ij2 c n (A, r) c n (/ 2 , r) (G(^o) + G(z n , A o))) 



n=l 



1 / °° 

^ 2 = 6 m ^ Sjir S,rA - F M + 2 £ e" n2 



c n (f 2l r)c n (f 1 ,r) 



x (F(^ fc o) - F(z n>A o))) + m h Zn J2T Znj 1t (G{z h o) + G{z A o) 

oo \ 

+ Z^e^^^K/.rKlA.rKGM + Gfv))) , (10) 



n=l 



2 Notice that, in this scheme, the self energy diagrams for on-shell leptons vanish since they can be written as 
XXf) = (P~ m fi)J2{p)(P~ m f2) 7 however, the vertex diagrams a and b in Fig. [I] give non-zero contribution. In 
this case, the divergences can be eliminated by introducing a counter term with the relation V^ en = V®+V^ , 
where V^ en (V®) is the renormalized (bare) vertex and by using the gauge invariance: k^V^ en = 0. Here, fc M 
is the four momentum vector of the outgoing photon. 

3 We take into account only the internal r-lepton contribution since, we respect the Sher scenerio [35], results 
in the couplings £jy y (i, j = e, /i), are small compared to £,§ Ti (i = e, fi, r), due to the possible proportionality 
of them to the masses of leptons under consideration in the vertices. Here, we use the dimensionful coupling 



%N,ij wrtn tne definition ^N,ij = & y wnere N denotes the word "neutral" 
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for fx (fz) = t; fi(fi or e; e). Here C\ = F % 2 ™a h , At (A 2 ) is the left (right) chiral amplitude, 
Zg — , z n g = 2 Ij t m\i i Qr is the charge of r lepton and the parameters c n , c' read 

771. ^ ' 171, g ~\~ yflj f -K ) 

c„(/,r) = cos[ "^ + ^) ], 

<U,r) = cos[rLfe|±M]. (11) 

In eq. f TTOT) the functions Fx(w), F 2 {w) are given by 

w(3-4w + «) 2 + 2lnw) 



(3 



w (2 + 3 w — 6 w 2 + w 3 + 6 w In w) 
G(u>) = - - ■ — . (12) 

(-1 + w) 4 

Now, we present the amplitudes A\ and A 2 appearing in the decay width T of the radiative 
LFV decays fi — > / 2 7 (see eq. ([9])) in the case of two extra dimensions (see Appendix section 
for details) 

x c 2 („, s) (/ 2 , t) (F(z {njS)ih o) - F(2; (ri)S)jA o))J + m fl t=§* rf2 £^ r/i (G(^o) + G(z A o) 



+ 4 E e- (n +s )<J /2i? c 2(n , s) (/i,r)c 2(n , s) (/ 2 ,r) (G(z (n , a))ft o) + G?(z (n>a) ^o)) y 

1 / 00 

^ 2 = ^J^[ Qm ^N, h r^r h (FM-FM+4Y, e- (n2+s2)a2/2R2 c' 2M (h,r) 
4o m T \ n )S 

x c 2 (n, s ) (A, r) (F(z (riiS)A o) - F(z (riiS) ^o))J + m /x £§ j2T £§* flT (G(z h o) + G(z A o) 

oo \ 

+ 4 £ e-(« 2 +^ 2 / 2 * 2 c' 2(niS) (/ 2 ,r)c 2(njS) (f^r) (G(z Mh o) + G(* (n> . )|A o))) , (13) 

n t s / 

2 

where the parameter Z( n ,s),s is defined as, 2( n ,s),s = m ^+ ra 2 7^2+52/^2 ■ I* 1 ec L- (O the summation 
would be done over n, s — 0, 1, 2..., except n = s = 0. Furthermore the parameters c 2 („ )S ) (/, r) 

and c 2(n, s ) (/» read 

4 M (/,r) = CM [ "fc/*+»ri+'(*/* + *«) 1 , (14) 

3 Discussion 

The radiative LFV decays /1 — > / 2 7 exist at least in the one loop level in the 2HDM where 
the tree level FCNC interactions are permitted. The Yukawa couplings £jv^-, i, j = e,fi,r, 
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which are the free parameters of the model, play an essential role on the physical parameters 
of these decays. Since we follow the idea that the hierarchy of lepton masses are due to 
the lepton Gaussian profiles in the extra dimensions, there appear exponential suppression 
factors, originated from the different locations of various flavors and their left and right handed 
parts of lepton fields, in the Yukawa part of the lagrangian, after the integration of the extra 
dimension(s) (see the eq. (JTj) (eq. (fT7|) ) for n = (n, s = 0)). We take the Yukawa couplings in 
four dimensions as the combination of these new factors and the higher dimensional one (see 
eq. (jHJ) and (118]) ) and consider that the couplings £jv,y; h3 = e '/ u are smaller compared to 
^fj ri i = e, /i, r since latter ones contain heavy flavor. Furthermore, we assume that, in four 
dimensions, the couplings is symmetric with respect to the indices i and j. 

Our analysis is devoted to the prediction of the effects of the extra dimensions on the LFV 
radiative decays. Here we choose the appropriate numerical values for the Yukawa couplings, 
by respecting the current experimental measurements of these decays (see Introduction section) 
and the muon anomalous magnetic moment (see [37j and references therein). Notice that, for 
the Yukawa coupling £§ TT , we use the numerical value which is greater than the upper limit 
of £jv t „. We also study the effects of the parameter p = a/R, where a is the Gaussian width 
of the fermions (see [18] for details). For the compactification scale 1/R, there are numerous 
constraints for a single extra dimension in the split fermion scenario. The direct limits from 
searching for KK gauge bosons imply 1/R > 800 GeV. The precision electro weak bounds 
on higher dimensional operators generated by KK exchange place a far more stringent limit 
1/R > 3.0 TeV [38]. In [22], the lower bounds for the scale 1/R have been obtained as 
1/R > 1.0 TeV from B — > (pK s , 1/R> 500 GeV from B -> ip K s and 1/R > 800 GeV from 
the upper limit of the BR, BR (B s — > p + p~) < 2.6 x 10~ 6 . We make our analysis by choosing 
an appropriate range for the compactification scale 1/R, by respecting these limits in the case 
of a single extra dimension. For two extra dimensions we used the same broad range for 1/R. 

In the model we use there are various free parameters, the new Yukawa couplings, the masses 
of new Higgs bosons and the ones coming from the split fermion scenario, namely, the com- 
pactification scale and the possible locations of fermions in the extra dimensions. It is obvious 
that the predictions are sensitive to those parameters and it is not easy to decide whether the 
enhancement comes from the broad region of one parameter set belonging to the 2HDM part, 
or the one belonging to the split fermion scenario. However, with the possible forthcoming 
experimental results of the processes which are more sensitive to the new parameters coming 
from the new Higgs doublet, the more stringent restrictions can be obtained. This would lead 
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to more accurate analysis of the effects due to the split fermion scenario. In our case, we try 
to estimate the sensitivity of the BRs to the split fermion scenario by fixing the other param- 
eters and we expect that the more accurate discussion of this scenario can be reached with 
forthcoming experiments which reduces the sensitivities at several orders. 

In the present work, we take split leptons in a single and two extra dimensions and use a 
possible set of locations to calculate the extra dimension contributions. We make the analysis 
in one and two extra dimensions. In the case of a single extra dimension (two extra dimensions) 
we use the estimated location of the leptons given in eq. ([3]) (eq. ( |T9l) ) to calculate the lepton- 
lepton-Higgs scalar KK mode vertices. For two extra dimensions, first, we take that the leptons 
are restricted to the fifth extra dimension, with non-zero Gaussian profiles. This is the case that 
the enhancement in the BRs of the present decays is relatively large. The reason beyond the 
enhancement is the well known KK mode abundance of Higgs fields. Finally, we assume that 
the leptons have non-zero Gaussian profiles also in the sixth dimension and using a possible 
set of locations in the fifth and sixth extra dimensions (see eq. (fT§j) ). we calculated BRs of the 
decays under consideration. In this case the additional the exponential factor appearing in the 
second summation further suppresses the BRs. 

In Fig. [2] ([3] ; BJ , we plot the BR of the decay \i — > (r — > ; r — > /jfy) with respect to the 
compactification scale 1/R for p = 0.01, m h o = 100 GeV, m^o = 200 GeV and the real couplings 
^ = lOGeV ^ >Te = 0.001 GeV (C§ ;TT = 100 GeV, £§ >Te = I GeV ; £* TT = 100 GeV, 
£n,th = 10 GeVjj. Here the solid (dashed, small dashed, dotted) line represents the BR without 
extra dimension (with a single extra dimension, with two extra dimensions where the leptons 
have non-zero Gaussian profiles in the fifth extra dimension, with two extra dimensions where 
the leptons have non-zero Gaussian profiles in both extra dimensions) . It is observed that BR is 
weakly sensitive to the parameter 1/R for the 1/R> 500 GeV for a single extra dimension. The 
enhancement of the BR is relatively larger for two extra dimensions due to the Higgs scalar KK 
mode abundances and, in this case, the sensitivity is weak for 1/R> 1.0 TeV. However, these 
contributions do not increase extremely due to the suppression exponential factor appearing in 
the summations. Furthermore, the numerical values of BRs are slightly greater in the case that 
the leptons have non-zero Gaussian profiles only in the fifth extra dimension. 

Now, we would like to present the amount of the enhancements in the BRs of the decays 
we study, by taking the existing bounds of the compactification scale 1/R into account and to 



4 For r — > e"f we take the numerical value of the coupling Te , f]| Te = 1 GeV. Here we try to reach the 
new experimental result of the BR of this decay (see j^j)- With the more sensitive future measurements of the 
BRs of these decays these couplings would be fixed more accurately. 
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discuss the possibility of observations of these additional contributions. For p — > (r — > e'j ; 
r — > /i7) decay, the enhancement in the BR for a single extra dimension is at the order of 0.7% 
(1.0% ; 1.3%) for 1/R ~ 800 GeV, compared to the case where there is no extra dimension. 
For the greater value of the scale 1/R, 1/R ~ 3.0 TeV, the enhancement reads 0.01% (0.03% ; 
0.05%). These numbers show that the enhancement in the BR of the decay r — > p^ is larger 
compared to the BRs of the others and this decay is probably the better candidate among the 
present LFV decays to detect the effects of the extra dimensions for a single extra dimension. 
Notice that, for 1/R ~ 3.0 TeV, the enhancements in the BRs of these decays are weak and 
difficult to observe. 

For two extra dimensions, where the leptons have non-zero Gaussian profiles in the fifth 
extra dimension, the enhancement in the BR of the decay p — > c-f (r — > e-y ; r — > /ry) is at 
the order of the magnitude of 3.6% (9.5% ; 10.6%) for 1/R ~ 800 GeV, compared to the case 
that there is no extra dimension. For 1/R ~ 3.0 TeV, the enhancement reads 0.05% (0.5% ; 
0.6%) ^. This shows that the additional contributions cause to increase the BR at the order 
of ~ 10% for t — > e7 and r — > p'y decays, in the case of 1/R ~ 800 GeV. Therefore, for 
two extra dimensions, r — > e^ and r —>■ p'j decays can ensure valuable information on the 
effects of extra dimensions in the split fermion scenario. Notice that, the lower bound of the 
compactification scale probably will be different in the case of two extra dimensions, compared 
to the one existing for a single extra dimension. However, in our calculations, we used the same 
broad range of the scale 1/R for a single and two extra dimensions. 

Fig. [5] (0 ; [7]) is devoted to the parameter p dependence of the BR of the decay p — + ej 
(r — > e7 ; r — > pj) for 1/R = 500 GeV, m h o = 100 GeV, ttiao = 200 GeV and the real couplings 
e| jTM = 10 GeV, el, re = 0.001 GeV (£$ >TT = 100 GeV, £# )Te = lGeV ; £f )TT = 100 GeV, 
Cn,tu, = 10 GeV). Here the solid (dashed, small dashed, dotted) line represents the BR without 
extra dimension (with a single extra dimension, with two extra dimensions where the leptons 
have non-zero Gaussian profiles in the fifth extra dimension, with two extra dimensions where 
the leptons have non-zero Gaussian profiles in both extra dimensions). The BR of the decay 
/i — > e7 (r — > e-y ; r — > fij) changes 0.8% (1.0% ; 0.8%) in the given interval of the parameter 
p and for 1/R = 500 GeV, in the case of a single extra dimension. This shows that the BRs 
of the present decays are weakly sensitive to p, which is the main parameter controlling the 
Gaussian widths and the possible locations of the fermions in the extra dimensions. For two 



5 For two extra dimensions where the leptons have non-zero Gaussian profiles in both extra dimensions, the 
enhancement in the BR is at the order of the magnitude of 3.5% (8.6% ; 9.7%) for 1/R ~ 800 GeV, for /1 — > 
(t — ► cy f t — ► ^7) decay, compared to the case that there is no extra dimension. For 1/R ~ 3.0 TeV, the 
enhancement reads 0.05% (0.4% ; 0.5%). 
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extra dimensions, the sensitivity of the BRs to the parameter p increases. For the decay p — > cy 
(r — > ey ; r — > /ry) the BR changes 3.6% (19.0% ; 15.0%) in the given interval of the parameter 
p and for 1/R = 500 GeV, in the case of two extra dimensions where the leptons have non-zero 
Gaussian profiles in the fifth extra dimension. For two extra dimensions where the leptons have 
non-zero Gaussian profiles in both extra dimensions, the BR of the decay p — > ey (r — > ey ; 
r — > p'-f) changes 5.1% (24.0% ; 18.0%) in the given interval of the parameter p. It is shown 
that the sensitivities of the BRs of studied LFV decays are relatively greater for two extra 
dimensions and this sensitivity increases for the case where the leptons have non-zero Gaussian 
profiles in both extra dimensions. In addition to this, the BRs of tau decays are sensitive to 
the parameter p and this sensitivity increases for p < 0.05, in the case of two extra dimensions. 

At this stage we would like to discus the possibilities of detecting the additional contri- 
butions due to the extra dimensions with the present and possible forthcoming experimental 
measurements. The experimental work for the lepton flavor violating decays has been done 
since the discoveries of heavy leptons. A new experiment, to search for the lepton flavor vi- 
olating decay p — >• ey [39j at PSI has been described and the aim of the experiment was to 
reach to a sensitivity of BR = 10~ 14 , improved by three order of magnitudes with respect to 
previous searches. At present the experiment (PSI-R-99-05 Experiment) to search the p — > ey 
decay is still running in the MEG [ID]. For the r — > /17 decay, recently, an upper limit of 
BR = 9.0 (6.8) 10 -9 at 90% CL has been obtained [H] ([12]) and this result is an improvement 
almost by one order of magnitudes with respect to previous one. The future measurement of the 
radiative p — > ey decay with the sensitivity of BR = 10~ 14 , hopefully, would make it possible 
to detect possible the additional contributions, especially, coming from two extra dimensions, 
even for the scale 1/R ~ 3.0 TeV. On the other hand, for the decay r — > py, the enhancement 
in the BR in the case of two extra dimensions is at the order of the magnitude of ~ 0.6% 
for 1/R ~ 3.0 TeV, and the one order improved experimental value of the BR would ensure a 
possible detection of the extra dimension effects even for large values of the compactification 
scale 1/R. 

As a summary, the BR is weakly sensitive to the parameter 1/R for 1/R> 500 GeV for a 
single extra dimension, however, this sensitivity increases for two extra dimensions. The expo- 
nential suppression factor appearing in the summations reduces the extreme enhancement due 
to the Higgs KK mode abundance. Furthermore, the BR is weakly sensitive to the parameter 
p especially for a single extra dimension case. In the two extra dimensions, this sensitivity is 
slightly larger compared to the one in the single extra dimension. With the help of the forth- 
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coming most accurate experimental measurements of the radiative LFV decays, especially the 
r — > /ry decay, the valuable information can be obtained to detect the effects due to the extra 
dimensions in the case of split fermion scenario. 
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5 Appendix 



In the two extra dimensions, after the compactification on the orbifold (S 1 x S 1 )/Z2, the new 
Higgs field 02 is expanded as 

1 ( oo ~j 

<h(x,V,z) = ^^°\x) + 2J2<P ( 2 ' S \x)cos(ny/R + sz/R)j . (15) 
The KK modes of charged (neutral CP even, neutral CP odd) Higgs fields existing in the new 



Higgs doublet have the masses Jm H± + + m? s , (Jm^ + + m 



m 2 A0 + + m? s ) where m n = n/R and m s = s/R are the masses of n'th and s'th level KK 
particles. If the leptons are restricted only the fifth dimension, the vertex factor Vlr(rl)%3 * s 
the same as the one in eq. (CO). In the case that the leptons are accessible two both dimensions 
with Gaussian profiles 

l lL = iVe-( ( ™ )2+(z -^ )2 ) /2CT2 / lL , 
E jR = Ne'( {m)2+{z - z ^ 2 ) /2(72 E 3R , (16) 

and the normalization factor N = 1 y 2 , the integration of the part of the Lagrangian 
fiL(R) S( n ' s \x) cos(ny / R + sz/R) fj R (L), with S = h°, A , over the fifth and sixth extra dimen- 
sions results in the vertex factor 

T/(n,«) _ ~{n 2 +s 2 )a 2 /4R 2 ~ ( (ViL (R) -VjR(L)) 2 +(ZiL (H) -z jR (£)) 2 ) /4ff 2 

LR (RL) ij ~ e e 

< n (y iL{R) + y jR{L) ) + s (z iL(R) + z jR(L) ) 

1 2R J ' 1 J 

Similar to a single extra dimension case, we define the Yukawa couplings in four dimension as 

*§ ((*?)*) = V lr1r L) y tf« /27TJ2 , (18) 

where V^Ln ■• = e _( ^Mfl)-%B(i)) 2 -*iMfl) V 4 ^ 2 . 

Lilt yri_L) %y 



Now we present a possible positions of 
dimensions, by using the observed masses 



eft handed and right handed leptons in the two extra 
1 

j. With the assumption that the lepton mass matrix 
is diagonal, one of the possible set of locations for the Gaussian peaks of the lepton fields in 
the two extra dimensions can be obtained as 



PL = V2 



a 



I (8.417,8.417) 

(1.0,1.0) ] , P ei = V2a 
\ (0.0,0.0) 



/ (4.7913,4.7913) 

(3.7913,3.7913) j . (19) 
V (-2.2272,-2.2272) 



3 The calculation is similar to the one presented in [18j which is done for a single extra dimension. 
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where the numbers in the parenthesis denote the y and z coordinates of the location of the 
Gaussian peaks of lepton flavors in the extra dimensions. Here we choose the same numbers 
for the y and z locations of the Gaussian peaks. 
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Figure 1: One loop diagrams contribute to l\ — > Z27 decay due to the zero mode (KK mode) 
neutral Higgs bosons h° and A (h 0n and A 0n ) in the 2HDM, for a single extra dimension. 
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Figure 2: BR(fi — > cy) with respect to the scale 1/R for p = 0.01, m^o = 100 GeV, m^o = 
200 GeV and the real couplings £f ;T/1 = lOGeF, ^ re = 0.001 GeV. Here the solid (dashed, 
small dashed, dotted) line represents the BR without extra dimension (with a single extra 
dimension, with two extra dimensions where the leptons have non-zero Gaussian profiles in the 
fifth extra dimension, with two extra dimensions where the leptons have non-zero Gaussian 
profiles in both extra dimensions). 
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Figure 3: BR(r — > ej) with respect to the scale 1/R for p = 0.01, m h o = 100 GeV, m A o = 
200 GeV and the real couplings £§ TT = 100 GeV, £§ Te = I GeV. Here the solid (dashed, small 
dashed, dotted) line represents the BR without extra dimension (with a single extra dimension, 
with two extra dimensions where the leptons have non-zero Gaussian profiles in the fifth extra 
dimension, with two extra dimensions where the leptons have non-zero Gaussian profiles in 
both extra dimensions). 
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Figure 4: BR(r — > /ry) with respect to the scale l/R for p = 0.01, m h o = 100 GeV, m^o = 
200 GeV and the real couplings £§ TT = 100 GeV, £§ T „ = 10 GeV. Here the solid (dashed, small 
dashed, dotted) line represents the BR without extra dimension (with a single extra dimension, 
with two extra dimensions where the leptons have non-zero Gaussian profiles in the fifth extra 
dimension, with two extra dimensions where the leptons have non-zero Gaussian profiles in 
both extra dimensions). 
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Figure 5: The same as Fig. [2] but with respect to parameter p and for 1/R = 500 GeV. 
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Figure 6: The same as Fig. [3] but with respect to parameter p and for 1/R = 500 GeV. 
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Figure 7: The same as Fig. H] but with respect to parameter p and for 1/R = 500 GeV . 
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